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Ab initio techniques are used to calculate the effective work function (Wefr) of a 
TiN/Hf02/Si02/Si stack representing a metal-oxide-semiconductor (MOS) transistor 
gate taking into account first order many body effects. The required band offsets were 
calculated at each interface varying its composition. Finally the transitivity of LDA 
calculated bulk band lineups were used and completed by MBPT bulk corrections 
for the terminating materials (Si and TiN) of the MOS stack. With these corrections 
the ab initio calculations predict a Wefr of a TiN metal gate on Hf02 to be close to 
5.0 eV. 
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Reducing the scale of the Metal-Oxide-Semiconductor (MOS) transistor has led the semi- 
conductor industry to major changes in the MOS stack. First, the Si02 oxide has been re- 
placed by a Hf02 layer on top of a Si02 thin film, while TiN is a leading candidate to replace 
poly-silicon in the gate for the next generation of CMOS transistor (22 nm and sub-22 nm). 
However, the control of the threshold voltage remains one of the major issues. The effective 
work function (Wcg) of a metal in a MOS structure is one of the key properties governing 
this threshold voltage. The aim of this paper is to construct a milestone for the comparison 
between an ideal system and available experimental data. 

During electrical capacitance-voltage (C-V) measurements Weff is usually estimated 
through: 

W,s = <Pms + X{Sl)+aSl). (1) 

Here xi^i) is the electronic affinity of silicon, C{Si) is the difference between the conduction 
band and the Fermi energy of doped silicon, and (pms is the difference between the metal 
Fermi energy and the semiconductor Fermi energy, (p^s is extracted from the flat-band 
voltage through modeling of the defect's charge distribution along the oxide^. 

The valence band offset (VBO) between the metal and the semiconductor (VBOms) is 
related to (p-^ns through VBO^s + Eg{Si) = (pms + CiSi), where Eg(Si)(l.l eV) is the electronic 
gap of silicon. Combining this relation with equation 1 and considering the well-accepted 
value of xiSi) (4.1 eV), the ab initio Wcs of a given metal can be evaluated through the 
equation 

^ab initio ^ ygQab initio ^ ^ ^^^g-y (2) 

Then one realizes that the ab initio quantity of primary importance is VBOms- 

Evaluating the VBOms from first principles calculations is a challenge. On the one hand 
the recent introduction of hafnium oxide requires the accurate treatment of the complex 
chemistry of the interfaces. On the other hand, in order to reach a sufficient precision on 
electronic energies (±0.1 eV), one needs to compute the eigenvalues in Many Body Pertur- 
bation Theory (MBPT). However, this is out of reach for the treatment of a complete MOS 
stack within the GW approximation (GWA, where G stands for the Green's function and W 
for the screened potential). Fortunately, the interface's potential drop can be well estimated 
within the Density Functional Theory (DFT)^, and coupled with the GW eigenvalues for the 
valence band states calculated only for each bulk material. These eigenvalues are calculated 
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in reference to the electrostatic potential average^ of each bulk. This reference potential 
is then aligned with the potential drop at the interface producing the band offset^. This 
technique is generally applicable for a sufficient material thickness (1 to 2 nm thickness for 
each slab) such that each material approaches a bulk state away from the interface. 

In this manner, the first order implementation of GWA, called GqWo, has been applied to 
correct the band alignments along Ill-V heteroj unctions. These corrections slightly improved 
the DFT results in the local density approximation (LDA), which were already in pretty 
good agreement with experiments^'^. Recently, this ab initio method has been applied 
to assess the band alignment at a junction between an oxide (Si02, Zr02 or Hf02) and a 
semiconductor (Si) resulting in a close agreement with experiment unlike the LDA results^'^. 
In this paper we use the same approach for each interface of a complete TiN / Hf02 / Si02 
/ Si stack. 

The potential drops at each of the interfaces have been computed within DFT/LDA using 
the atomic orbital SIESTA software^ whereas the GWA computations have been carried 
out with the plane- waves ABINIT software^. To be consistent, identical TrouUier-Martins 
pseudopotentials^ for each atomic species have been used for Abinit and Siesta and were 
generated using the Fhi98pp and Atom programs^'^^ in the local density approximation 
(LDA)^^. Siesta calculations have been performed using a polarized double zeta, DZP, basis 
with an energy shift of 50 meV and a Meshcutoff of at least 100 Hartrees. A simulated 
anneal has also been carried out within Siesta as an heuristic method to relax the forces for 
each stack: from lOOOK to lOOK in 1 picosecond, followed by a conjugated gradient method 
to reduce the forces to less than 0.07 eV/A. 

All GWA calculations were performed at the GqWo level. The quasiparticles energies 
were carefully converged with respect to the energy cut-offs, the number of k-points in the 
Brillouin zone, and the number of bands used to compute the dielectric function and the self- 
energy. Thus a numerical accuracy of ±0.1 eV on the quasiparticles energies was reached. 
The Plasmon-Pole model proposed by Godby and Needs^^ was used to describe the dynamic 
dependence of the "screening" function. 

Two types of Si/Si02 interfaces have each been simulated, one oxygen poor and one 
oxygen rich. These supercells contained a 22.2 A long slab of /3-cristobalite Si02 and a 24.3 
A long slab of diamond Si, both oriented with the (001) plane at the interface. The latter 
slab was compressed by 3% in the [100] and [010] directions. The simulated LDA value 
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of the VBO^^"^ is +3.0 eV (+2.8 eV) for the oxygen rich (poor) interface using the sign 
convention where VBO>0 when the Valence Band Maximum (VBM) of the material on the 
left is above the VBM of the material on the right. The GqWo corrections to the VBMs 
computed for each material are 5e'^'^~^^'^{Si) = —0.6 eV and 6e^^~^^^{Si02) = —1.9 
eV, which lead to a VBOf^ of +4.2 ± O.leV. These results are in good agreement with 
the results published by Shaltaf et al? (4.1 eV calculated within GWA starting from GGA 
electronic structures), and a little lower than experimental data ranging from 4.3 eV to 4.5 
eV estimated by X-ray photoemission spectroscopy^^'^^. 





axis perpendicular to the interface (bohrs) 
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FIG. 1: O rich Si02/Hf02 interface (a) and nitridized TiN/Hf02 interface (c): N in purple, 

Ti in turquoise, Hf in blue. Si in pink and O in red. Potential projected along the 
Hf02/Si02(b) and the TiN/Hf 02(d) superlattices (black line) and smoothed by a double 
convolution (red line). The black arrows indicate the shifts of the eigenenergies with 
respect to the averaged potential in the bulk. The blue arrow indicates the VBO^^"^. 



The Si02/Hf02 interface is the second interface of the complete stack (figure 1) considered 
here. The Hf02 slab was simulated in its monoclinic structure with P2i/ c symmetry, which 
is its stablest phase at ambient pressure and temperature^^'^^ with the (001) plane at the 
interface. The Si02 slab was again /3-cristobalite and oriented with the (001) plane at the 
interface. The Hf02 was stretched by 3.7% in the [100] direction and by 2.1% in the [010] 
direction to match the 5.23 x 5.23 square interface of the Si02, which agrees with the other 
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simulations in literature of Si02 in the stack^'^^. The supercell contained seven Hf02 layers 
and ten Si02 tetrahedron. At the interface after relaxation, 3 oxygens were each three times 
coordinated with both Si and Hf. The potential along the stack shown in figure 1 exhibits 
a large potential drop of 5.2 eV, which is expected since the electric permittivity and the 
average valence electronic density are much higher in hafnia than in silica. At the LDA level 
the VBOg^"^ is +0.4 eV, which is qualitatively wrong with respect to XPS measurements^^, 
i.e. the VBM of Si02 is measured below the VBM of Hf02. Some DFT calculations of the 
Hf02/Si02 interface at the LDA level have also been performed by Sharia et al. showing 
that the valence band offset is highly dependent on the interface stoichiometry^^. Using our 
GqWo correction for Si02 (—1.9 eV) and for Hf02 (—0.5 eV close to the value of Griining et 
al.^ of —0.4 eV) the LDA band offset is reversed producing a VBO^^ of —1.0 eV. This is 
qualitatively and quantitatively in good agreement with XPS measurements^'^ of 1.05 ± 0.1 
eV. 

The TiN/Hf02 interface is the final interface of the stack. Constructing this interface 
requires a specific procedure because of the large lattice mismatch. First, a slab of 6 layers 
of monoclinic Hf02 surrounded by vacuum oriented in the direction perpendicular to the 
(001) plane containing three unit cells in the [100] direction and one unit cell in the [010] 
direction was generated. During molecular dynamic simulations several TiN layers were 
successively added to the surface. We found that the positions of the interfacial atoms from 
TiN were influenced by the underlying Hf02 slab arranging themselves in almost the same 
structure as the oxide (at ± 0.3 A of the atom sites). Thanks to this observation when 
we constructed our TiN/Hf02 stack model the Ti atoms at the interface were placed in 
the sites of the hafnium and the nitrogens in the oxygen sites. The stack was finalized by 
adding a face-centered-cubic TiN slab oriented with the (111) surface plane on top of the 
interfacial TiN layer, which minimizes the mismatch between the slabs. Indeed, it has been 
pointed out experimentally that the (111) growth direction lowers the film strain^^'^^. In our 
case the TiN is 1.8% stretched out in the [110] direction and 0.5% compressed in the [112] 
direction, and contained 6 TiN layers expanded by 1% in the [111] direction in order to keep 
the average electronic density of the TiN bulk. Two types of interfaces were computed: one 
"nitridized" with oxygen sites of Hf02 substituted by N at the interface (figure 1), the second 
"oxidized" , with six O replacing the six N at the interface. The VBOg^^ obtained for these 
two interfaces are —2.4 eV (figure 1) and —2.5 eV respectively. Previous DFT calculations 
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TABLE I: Valence band (VBO) and conduction band (GEO) offsets with the resulting 
effective work function. The sign of the band offset is positive when the VBM of the 
material on the left is above the VBM of the material on the right. 



interfaces Si/Si02 


Si02/Hf02 


Hf02/TiN 


Si/TiN 


Weff(TiN) 


rich poor 




nitridized oxidized 






VBO(LDA) +3.0 +2.8 


+0.4 


-2.4 


-2.5 


[+0.7, +1.0] 


[5.9, 6.2] 


VBO(GoWo) +4.3 +4.1 


-1.0 


-3.3 


-3.4 


[-0.3, 0.0] 


[4.9, 5.2] 


CBO(GoWo) -3.2 -3.4 


+1.6 


+2.7 


+2.6 


[+0.8, +1.1] 





have been done on TiN/Hf02 stacks but cannot be compared with ours since they were 
computed in the generalized gradient approximation^*^. The quasiparticle correction to the 
IDA Fermi energy of TiN was calculated to be Se^^-^^^{TiN) = +0.4 eV. 

The valence band alignments along the whole stack at the LDA and GqWo levels where, 

VBO„,s = VBOi + VBO2 + VBO3, (3) 

are depicted in table I. Also included in table I are the conduction band alignments cal- 
culated by using our GqWo band gaps for Si (1.1 eV), Si02 (8.6 eV), and Hf02 (6.0 eV). 
The VBOji^^^ ranges between —1.0 eV and —0.7 eV depending on the chemistry of the in- 
terfaces which would correspond to a work function of between 6.2 eV and 5.9 eV. This is 
in qualitative disagreement with the experimental Wcfr which is usually equal to or below 5 
eV. It is worthy of note that in order to calculate the VBO between the two materials on 
the extremities of the stack one only needs to compute the GW energies for the top and the 
bottom slabs because for the intermediate slabs the corrections vanish yielding, 

VBO^w = VBO^DA + Se^'^-'^^'^iSt) + 6e^'^-'^''^{TtN). (4) 

VBO^^ ranges from 0.0 eV to 0.3 eV depending on the chemistry of the interfaces, and 
one finally obtains an ab initio W^s evaluated in the MBPT framework ranging from 5.2 eV 
to 4.9 eV (i.e. P+ in jargon of the semiconductors industry). The Wefr for (111) oriented 
TiN is in agreement with those measured during cold growth process which mostly does not 
modify the TiN stoichiometry^'^^. We also point out that the GW corrections applied with 
a coherent methodology are mandatory in order to obtain qualitatively correct results for 
the evaluations of VBOs and W^s in a MOS stack. 
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We think that our results can be used to adjust the parameters of the model used to 
simulate the whole MOS device. Secondly, by comparing our 'ideal' model system (no 
defects and a perfect stoichiometry) with the dispersion of the experimental values (from 
4.2 to 4.9 eV) assessed by C-V measurements^ '^^ we think that our values can be used as 
a reference to study the impact of the fine chemistry in such a gate stack. Finally, in this 
paper we show that the Wefr is an accessible and reliable ab initio quantity and that this 
type of calculation is promising to help the understanding of the Wefr for nanodevices in 
microelectronics. 
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